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ABSTRACT 

We describe a new multi-colour proper motion survey for cool white dwarfs (CWDs). 
The observational database consists of ~ 300 digitally scanned Schmidt plates in 
ESO/SERC field 287. The entire survey procedure, from the raw Schmidt plate data 
to final white dwarf luminosity function (WDLF) is described, with special emphasis 
on completeness concerns. 

We obtain a sample of 58 WDs, for which we have follow up CCD photometry 
and spectroscopy of a representative sub-sample. Effective temperatures and luminosi- 
ties of our sample objects are determined by comparing photometry with the model 
atmosphere predictions of Bergeron, Saumon and Wesemael. Space densities are cal- 
culated using an adaptation of Schmidts (l/Vmax) method, from which a WDLF is 
constructed. Comparison of our observational LP with the models of both Wood and 
Garcia-Berro et al. indicate an age for the local Galactic Disc of lO^g Gyr. Impor- 
tantly, we find no evidence of incompleteness in our survey sample. Proper motion 
number counts imply the survey is complete, and the WD sample passes the {V/Vmax) 
completeness test. 
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1 INTRODUCTION 

A reliable age for the local Galactic Disc places a valuable 
constraint on both globular cluster ages and cosmological 
models. A number of independent methods of investigating 
this problem have been employed in the past (eg. Jimenez 
1998 and references therein), resulting in a broad consen- 
sus that the lower limit for the Disc age lies between 8 and 
12 Gyr. Potentially one of the most reliable means of esti- 
mating the Disc age is via cool white dwarf (CWD) stars. 
These estimates use the idea, first proposed by Schmidt 
(1959), that in a galaxy of finite age there will be a temper- 
ature beyond which the oldest, coolest white dwarfs (WDs) 
have not had time to cool. This predicted cut-off in the lu- 
minosity function (LF) of WDs, if satisfactorily observed, 
can then be used in conjunction with WD cooling models to 
derive the Disc age. 

CWDs are difficult to find, being both extremely faint 
and of similar colour to the numerous K and M-type dwarfs; 
and are almost exclusively discovered by means of their 
proper motion. The cut-off in the WDLF was observed 
(Liebert, Dahn & Monet 1988, hereafter LDM) after thor- 
ough follow up observations of CWD candidates drawn from 
the Luyten Half Second (LHS) Catalogue (Luyten 1979). 
Although at that time a Disc age of 9.3 ± 2 Gyr was de- 



rived from this sample (Winget et al. 1987), further ob- 
servations and improvements in model atmospheres (Berg- 
eron, Ruiz & Leggett 1997, hereafter BRL) and theoreti- 
cal LPs (Wood 1992, 1995) has prompted a recent rede- 
termination of the Disc age for the same sample (Leggett, 
Ruiz & Bergeron 1998, hereafter LRB), yielding a value of 
8 ± 1.5 Gyr. While the existence of the cut-off in the LDM 
WDLF has not been challenged by subsequent observational 
work, the details of its precise position and shape have. A 
sample of CWDs found using common proper motion bina- 
ries (CPMBs), again culled from the Luyten surveys, suggest 
that there are ~ 5 times more very faint CWDs than found 
by LDM (Oswah et al. 1996, hereafter OSWH). A Disc age 
of 9.5 ± 1 Gyr was found using this sample, and the factor 
of ~ 5 increase in the faintest WDs has been confirmed by 
an independent search for CWDs in the south (Ruiz and 
Takamiya 1995). 

Until now, the proper motion catalogues used to ex- 
tract samples of CWDs have been produced by 'blink' com- 
parison. While these surveys have clearly been successful 
in picking up individual stars of low luminosity and high 
proper motion, the use of such a subjective survey tech- 
nique raises worries concerning completeness. The advent of 
high precision micro-densitometers such as SuperCOSMOS 
(Hambly 1998 and references therein) allow proper motions 
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and magnitudes to be calculated objectively using a series 
of plates in the same field. Hambly et al. (1997) have re- 
cently discovered in Taurus perhaps the coolest known WD 
using just such a procedure. This object, WD 0346+246, 
should certainly appear in the LHS Catalogue, and would 
then presumably have been included in the LDM sample 
since it fulfills the given criteria. This object was discovered 
serendipitously from work in a particular 6x6° Schmidt 
field, and therefore casts further doubt on the completeness 
of the Luyten catalogue. 

In this work we exploit both the extensive collection of 
over 300 plates in ESO/SERC field 287 and the power of 
the SuperCOSMOS measuring machine to produce a deep, 
complete multi-colour proper motion survey from which a 
sample of CWDs has been extracted. Much attention was 
given to the critical question of completeness, with specific 
reference to the choice of survey limits and the sources of 
potential contaminants. The paper is organised as follows: 
the plate database and reduction of the digitally scanned 
data is described in Section ^; the choice of survey limits 
to combat sample contamination and the method of WD 
sample selection are discussed in Section ^ in Section ^ the 
important question of the high proper motion limit is ad- 
dressed by independent tests; some follow up observations 
of our WD candidates are presented in Section ^ stellar pa- 
rameters are calculated for our sample in Section ^ which 
also includes a discussion of potential contaminant popula- 
tions; the WDLF from our sample is presented in Section ^ 
and the derived Disc ages discussed in Section H. 
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Figure 1. Number of plates per year in ESO/SERC Field 287 in 
the Bj (top) and Rp (bottom) passbands 



Year 


Bj 


R 


Year 


B,T 


R 


Year 


Bj 


R 


1977 


1 




1984 


4 


2 


1991 


1 


1 


1978 


1 




1985 


2 


1 


1992 


1 


1 


1979 






1986 


1 


1 


1993 


1 


1 


1980 




5 


1987 


1 


1 


1994 


1 


1 


1981 






1988 


1 


1 


1995 


1 


1 


1982 






1989 


1 


1 


1996 


1 


1 


1983 


1 


1 


1990 


1 


1 









2 SCHMIDT PLATE DATA REDUCTION 

In the course of a long term quasar variability study (eg. 
Veron & Hawkins 1995) over 300 Schmidt plates have been 
taken in ESO/SERC field 287. Figure | shows the distri- 
bution over time of the plate collection in this field in the 
two principal passbands available, Bj and Rf . Additional 
plate material exists in the U (19 plates), V (11 plates) and 
I (40 plates) passbands. All these plates have been digitally 
scanned by either the SuperCOSMOS plate measuring ma- 
chine or its predecessor COSMOS (MacGillivray and Stobie, 
1984). As described in Appendix a CWD survey utilis- 
ing the reduced proper motion (RPM) population discrim- 
ination technique requires photometry in at least two pass- 
bands and proper motion data. Although for most stellar 
types a survey of this sort will be proper motion limited, 
for extremely intrinsically faint objects such as cool degen- 
erates the photometric survey limits also become important. 
A principal concern in this work is therefore to maximise the 
survey depth. 

The technique of stacking digitised Schmidt plates has 
been in existence for some time (eg. Hawkins 1991; Kemp & 
Meaburn 1993; Schwartzenberg et al. 1996), and has recently 
been thoroughly investigated with SuperCOSMOS (Knox et 
al. 1998). The database presented in Figure |l]is an obvious 
candidate for stacking in a proper motion survey, since at 
least 4 plates are available in most years and the use of 4 
plate stacks rather than single plate data will yield deeper 
photometric survey limits with no reduction in the proper 
motion time baseline. The data has therefore been grouped 
into 4 plate stacks according to Table W 



Table 1. The number of Bj and R 4 plate stacks per annual 
epoch. 



In order to track objects through epochs and calculate 
relative proper motions, all the stacks have been shifted to a 
common co-ordinate system corresponding to a measure of 
the high quality Bj plate J3376. This is achieved by a global 
transformation followed by local transformations (transla- 
tion, rotation and scale), performed by splitting the stack 
area into a grid of 16x16 small areas (Hawkins 1986). Every 
object in each area is used to define the local transformation 
(since fainter objects are more numerous, they essentially 
define the astrometric coordinate system). Objects are then 
paired up between epochs. The software used in this work 
operates using a fixed 'box size' (13arcecs) in which it looks 
for a pair, thus in the first instance a high proper motion 
limit is imposed on the survey - this issue will be addressed 
in detail later. The pairing procedure yields a list of x and 
y coordinates for each object, one set for each stack the ob- 
ject is found on. Calculating proper motions is then simply a 
matter of performing a linear regression fit to each object's x 
and y coordinates as a function of time. However, erroneous 
pairing inevitably occurs between stacks, and we therefore 
wish to perform some form of bad point rejection to reduce 
contamination by spurious proper motions. In order to reject 
deviant points (and calculate parameters such as cr^ and x^) 
an estimate of the error associated with each measure of po- 
sition is required. We assume this error is simply a function 
of magnitude and that it will vary from stack to stack, but 
not across the survey area. This error is calculated using the 
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Figure 2. An example of a linear regression fit to a set of object 
positions. The vertical axis is deviation of measured position on a 
given stack from the mean object position, with year this century 
on the horizontal axis 



deviation of an object's position on a particular stack from 
the mean position over the 20 stacks used, and is determined 
over 10 magnitude bins. A Sa iterative rejection procedure 
is implemented to reject spurious pairings or high proper 
motion objects which are not reflecting the true positional 
errors sought. The calculated errors are much as one might 
expect: decreasing for brighter objects until factors such as 
saturation and blended images makes positional measures 
more uncertain. 

A straight hne fit can now be applied to the x and y 
data for each object, the gradient of which is taken to be the 
measured proper motion, and jiy respectively. An exam- 
ple is shown in Figure |2|, the points showing the deviation at 
each epoch from the average object position with error bars 
calculated as above. Deviant points arising from spurious 
pairings often lie far from the other data and will give rise to 
spurious high proper motion detections if not removed. We 
therefore iteratively remove points lying Scr from the fitted 
line. This can occasionally lead to further problems if there 
are several bad points associated with the object, and the 
result of several iterations can be a larger spurious motion 
detection. This source of contamination is generally elimi- 
nated by insisting sample objects are detected on virtually 
every stack. The validity of the positional error estimation 
scheme described above has been verified by confirming that 
scatter plots of log reduced a-s a function of magnitude 
cluster around zero for all magnitudes in all regions of the 
survey area. 

Instrumental magnitudes are calculated for every 
object detected on each stack in the standard COS- 
MOS/SuperCOSMOS fashion (Beard et al. 1990 and ref- 
erences therein). Briefiy, an object detection is defined by a 
given number of interconnected pixels with intensity above 
a given threshold (eg. 8 interconnected pixels with intensity 
above a 2.5a sky noise threshold for SuperCOSMOS data). 
An object's instrumental magnitude is then calculated as the 
log of the sum of the intensity above background across the 
object area. This quantity varies monotonically with true 
magnitude, and is therefore suitable for use in construct- 



Figure 3. Number counts as a function of apparent R magnitude 



ing calibration curves using a CCD sequence. A sequence 
of ~ 200 stars with CCD magnitudes measured in a vari- 
ety of passbands exists in field 287 (Hawkins et al. 1998) , 
yielding U, B, V, R and I photometry to a typical accuracy 
of 0.15 magnitudes (see Section 5.2). Significantly smaller 



errors are theoretically obtainable from photographic mate- 
rial, and the larger uncertainties we find appear to be caused 
by systematic deviations of sequence objects from the cal- 
ibration curve. This is not a colour or field effect, and is 
probably caused by differences in detection media. 



3 SURVEY LIMITS AND SAMPLE 
SELECTION 

The 'catalogue' resulting from the implementation of the 
procedure described in the previous section consists of astro- 
metric and photometric measures for over 200,000 objects. 
Criteria for inclusion in this preliminary sample is merely 
detection in both Bj and R passbands (since these are re- 
quired for construction of the reduced proper motion dia- 
gram (RPMD)) and a measure of proper motion in both 
these passbands. It is from these objects that an uncontam- 
inated proper motion sample is to be drawn; and we require 
well defined universal survey limits so that space densities 
can be calculated from the final survey sample. 

Number count plots from this survey data increase lin- 
early with increasing magnitude, as shown for the R data in 
Figure ^ before dropping precipitously. This cut-off is at- 
tributed to the survey detection limit, and the position of 
the turnover is used to determine photometric survey limits. 
The limits used are 21.2 in R and 22.5 in B. 

The proper motion distribution for all objects in our 
survey area detected on at least 15 stacks in both B and 
R is shown in Figure ^ Low proper motions are generally 
an artifact of measuring machine error, thus the distribu- 
tion indicates a typical error in measured proper motions of 
~ lOmas/yr. Our criteria for choosing a survey proper mo- 
tion limit are elimination of contaminant spurious proper 
motions from the sample and , with this in mind, maximis- 
ing the size of the final proper motion sample extracted. The 
deterioration of positional accuracy with magnitude means 
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Figure 4. Proper motion frequency distribution 
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that fainter objects also have more uncertain proper motion 
determinations. The peak of the proper motion distribution 
therefore moves to higher proper motions with object sam- 
ples drawn from progressively fainter magnitude bins. For 
this reason, we consider the survey proper motion limit as a 
function of magnitude. 

Three independent methods of determining an appro- 
priate proper motion limit have been investigated for this 
dataset: 

(1) analysis of proper motion error distribution 

(2) cumulative proper motion number counts 

(3) RPMD inspection 

all of which are described here. 

The characteristics of the 'noise' in the proper mo- 
tion distribution can be analysed by assuming all objects 
are 'zero proper motion objects' and the calculated proper 
motions arise purely from random measurement error. The 
'true' values of and measured by linear regression are 
therefore zero, and the random measurement errors give rise 
to a normal error distribution in and about zero with 
common a. We are of course interested in the total proper 
motion, which follows a Rayleigh distribution of the form 

P(p) = 4e-^. (1) 

This survey utilises two independent measures of proper mo- 
tion in the same field, one from the stacks of B plates and 
the other from the R stacks. A useful means of reducing the 
final proper motion survey limit will be to compare these 
two measures and reject inconsistent motions. Comparison 
of independent measures of proper motion needs to be in- 
corporated into this analysis if it is to be useful in predicting 
sample contamination later. The measured B and R motions 
and their associated error distributions are therefore char- 
acterised 



(2) 



and 

P{^Ir) = ^rAtre"'^'-^^ (3) 

respectively. The algorithm used to select sample objects on 



Figure 5. Objects will only be selected for the survey sample if 
their proper motions lie in the hatched region of the diagram. 



the basis of proper motions will use three sequential cuts 
to eliminate spurious motions. Firstly, the averaged proper 

ie: 



motion must exceed the proper motion limit /iHm , 

/ir + IJ.b 



(4) 



Secondly, the difi'erence in proper motion must not exceed a 
second survey parameter (A/i): 



|M6 - M'-l < (A/i), 



(5) 



and finally the difference in position angle must not exceed 
a third survey parameter (Af/i): 



^r\< (A0). 



(6) 



An object must satisfy all three of these criteria to be in- 
cluded in the sample. The desired outcome of this analysis is 
the ability to predict the expected contamination from spu- 
rious motions for a survey with parameters num, (A/x) and 
{A<j)) given the 'zero proper motion object' error distribu- 
tion described by equations |^ and ^ In order to do this we 
consider the combined probability distribution P{^i,)P{fir) 
on the fj.b - fJ.r plane. The survey criteria described above 
effectively limit the selected sample to a specific area on 
the Hi, - plane, shown in Figure |^. Given a realistic esti- 
mate of the (normalised) P{nb) and P{fJ,r) distributions an 
integration over the hatched region is Figure |^ yields an es- 
timate for the fraction of objects, X, belonging to the error 
distribution likely to contaminate the sample. The integral 

/ P{pb)PifJ-r)dflrdtJ,b + 

rMb + {AM) 

P{^ib)PifJ,r)dfirdflb (7) 

' J2 "'2^6 -(Am) 
where 



Ji = 
and 



(8) 
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Figure 6. Normalised B (dashed line and dots) and R (solid line 
and crosses) measured data(points) and fitted error distributions 
(lines). This example uses the first 15 data points to determine 
the error distribution. 



J2 



(9) 



is easily calculated numerically for arbitrary fium and (A/i). 
Since this analysis concerns zero proper motion objects 
whose spurious motion arises purely from machine measure- 
ment error, we may assume no preferred position angle, and 
the resulting fraction X will be cut by a further (2(A</>))/360 
((A(^) in degrees) by the position angle selection criterion, 
ie. 

" 2(A0)/36O- 

The calculation described above requires knowledge of the 
P(/ii,) and Pf/ii ) 'zero proper motion' error distributions. 
If equation m is to be used to calculate numbers of con- 
taminants it is important that and P{Hr) can be si- 
multaneously rescaled from normalised distributions. The B 
and R data are therefore paired so that P{jJ.b) and P{^J.r) 
contain the same objects and thus the same number of ob- 
jects. For the purposes of this analysis this proper motion 
data is thought of as consisting of two distinct distributions: 
the 'zero proper motion object' error distribution, on which 
a distribution of real proper motions is superposed. Any 
attempt to determine the error distribution from the mea- 
sured distribution must use only the low proper motion data 
where the random errors of interest dominate systematics 
introduced by real proper motions. In the example to be 
shown here the first 15 bins (Imas/yr bins) of the normalised 
proper motion distribution is taken to be representative of 
the error distribution. 

The first 15 data points in both B and R are fitted with 
an assumed error distribution like equation |^. These fits are 
shown along with the measured data in Figure |^, with the 
real data rising above the error fits representing real proper 
motions. Since the error distributions exclude real motions 
they must be renormalised before use in equation ^. The 
rescaling factor to be used in calculating numbers of con- 
taminants is taken to be the number of objects contained 
in the error distributions. This number will in general be 
different for the B and R data, so the average is used. Nu- 




Miim(mas/yr) 



Figure 7. The results of numerical integration predicting the 
number of contaminant objects in the survey sample as a function 
of survey proper motion limit. The lower of the lines used a (A/^) 
of 5 mas/yr, with (A/i) rising to 50 mas/yr for the top line ((A0) 
was set to 90 throughout). These calculations were based on error 
distributions determined from the first 15 proper motion bins of 
the measured fii, and /ir distributions. 



merical methods can now be used to predict the number of 
contaminant spurious proper motions for a range of iium, 
(A/i) and (A^) using equations |^ and |lol The result of a 
series of calculations is shown in Figure|7|. The predicted 
contamination falls rapidly below one object as /iHm exceeds 
~ 42mas/yr. Care must be taken however that the error dis- 
tributions are not overly sensitive to the number of points 
in the measured fib and fir distributions used in their cal- 
culation. The error distributions were therefore recalculated 
using the first 30 bins in the measured fib and fir distri- 
bution and the predicted contamination plot redrawn using 
these new distributions . The two calculations prediction of 
the fiiim at which the survey contamination drops to below 
1 object are consistent to within ~ 3mas/yr. Presented in 
Table ^ are the results of calculating the fium at which the 
predicted number of contaminant objects falls below one for 
a range of magnitude cuts, (A/x) and derived P{fib), P{fir) 
distributions. Table |^ indicates that the (Afi) survey param- 
eter has a small bearing on the fium where Ncontam ~ 1. 
The (A/i) criterion can therefore be relaxed substantially to 
ensure no real motions are rejected. A similar argument is 
applicable to the (A^) parameter, where the small advan- 
tage gained by tightening the (A^) criterion again becomes 
increasingly outweighed by the potential for rejection of real 
proper motions. To conclude comment on this analysis, any 
fiiim greater than ~ 45mas/yr should eliminate contamina- 
tion arising from normally distributed positional measure- 
ment errors, although thus far no account has been taken of 
spurious motions arising from other sources (eg. erroneous 
pairings) . 

The second means of investigating the proper motion 
limit is via proper motion number counts. Assuming the lo- 
cal Disc has constant stellar density and that Disc kinemat- 
ics give rise to an approximate inverse correlation between 
distance and proper motion, a plot of log cumulative num- 
ber count (from large to small fi) versus log fi should fol- 
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R magnitude No. of Hb^t^r Miim (mas/yr) for varying (Afi) 

cut data points used (A/i)=5 (A^)=10 (A/i)=30 (A/i)=50 



_R > 20.5 15 40 41 42 43 

25 41 43 44 44 

20.5 > i? > 20 15 34 35 36 36 

25 34 35 36 36 

20 > R > 17 15 30 31 32 32 

22 32 33 33 33 

i? < 17 15 31 32 33 33 

25 34 35 36 36 



Table 2. Calculated where predicted number of contaminants falls below 1 object. Column two represents the number of bins in 

the measured fii, and fir distributions used to derive the P{in,), P(fj,r) distributions. 
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Figure 8. Cumulative number counts plot showing deviation 
from simple model at ^ 50mas/yr. 




^ ^ 

-10 12 3 -10 12 3 

Bj-R Bj-R 

Figure 9. A RPMD of a sample with iiiim= 50mas/yr for all 
objects. 



low a straight line of gradient —3 (ie. log^ A'^ oc — 31og/x). 
In principle a proper motion limit for this survey could be 
obtained by determining the point at which our measured 
proper motions deviate from this relation due to the exis- 
tence of spurious motions. In Figure |^ the cumulative num- 
ber counts are plotted as a histogram. If the points between 
log/i = 1.9 and log/i = 2.5 are fit with a straight line the 
resulting gradient is (— 2.998±0.071), in excellent agreement 
with the idealised predicted slope of —3. The fit is shown as 
a dashed line. This finding compares favourably with a sim- 
ilar analysis of Luyten (1969, 1974, 1979) and Giclas (1971, 
1978) common proper motion binary stars, which on an anal- 
ogous plot describe a straight line of significantly shallower 
gradient that the expected —3 indicating increasing incom- 
pleteness with decreasing proper motion (Oswalt & Smith 
1995). Towards lower proper motions our data rises above 
the line, indicating the onset of contamination at /i 50 
mas/yr, lending credence to the findings of the previous 
analysis. 

The final means of assessing the effects of varying /iHm 
is the RPMD. The RPMD is, for all sample objects lacking 
follow up observations, the sole means of stellar population 
discrimination. For this reason it is worthwhile inspecting 
the RPMD of samples produced using various survey lim- 
its. The principal concern is that the white dwarf locus be 



a. 

O 
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Figure 10. A RPMD of a sample with fiiim= 80mas/yr for all 
objects. 



as distinct as possible at all colours while maximising the 
sample size. 

It can be argued that there is little purpose in low- 
ering fiiim for bright objects to anything near the level of 
potential contamination. Consider searching for an intrinsi- 
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cally faint star such as a reasonably cool white dwarf, with 
an absolute R magnitude of ~ 14, amongst survey objects 
with apparent R magnitudes as faint as ~ 19. Such an ob- 
ject, having a (conservative) tangential velocity of 40km/s, 
would have a proper motion of ~ 80mas/yr - well beyond 
the predicted contamination limit. This argument becomes 
even more forceful for both intrinsically fainter and appar- 
ently brighter objects, implying that a conservative nnm is 
desirable for all but the faintest survey objects. 

Figure ^ is a RPMD produced with a /iHm of 50 mas/yr 
at every magnitude. The main sequence and white dwarf loci 
are visible, as is a dense group of objects between the two at 
Hji ~ 19.5. This group of objects probably consists of two 
sub-groups. Firstly of course, the most likely sample con- 
taminants are faint objects creeping through just over the 
proper motion limit at this point. Secondly, one would ex- 
pect a large population of objects to lie on the RPMD where 
the n and R distributions are most populated , ie. at fiHm 
and at faint R (an object with jj, — 50 mas/yr and R=21 
has Hr = 19.5), and it should therefore not be automati- 
cally assumed that every object lying at this point in the 
RPMD is suspect. A cause for concern however, is the way 
this population bridges the gap between the main sequence 
and the white dwarf population, a property one would ex- 
pect from a contaminant locus rather than the detection of 
bona fide proper motions. If a more conservative limit of 
A'iim=60 mas/yr is adopted for objects with R > 20.5 the 
RPMD (Figure |ll|) looks more promising, with the 'con- 
taminant locus' all but gone, leaving only the expected mild 
confusion (Evans 1992) between populations at their faint 
extremity. The adoption of an extremely conservative fium 
of 80 mas/yr for all objects leads to an even more well de- 
fined RPMD (Figure UW, where the white dwarf population 
discrimination is almost without exception unambiguous. 

The survey parameters to be used to select the prelimi- 
nary sample (ie. a sample subject to further object-by-object 
scrutiny and potential rejection, with the possibility of fur- 
ther objects being included the sample which lie just outside 
the RPM cut on one or both RPMDs) have been chosen with 
reference to the findings in this section. 

Error analysis and number counts suggest a of 50 
mas/yr, a (A/x) of 50 mas/yr and a (Acp) of 90 degrees 
should essentially eliminate contamination arising from nor- 
mally distributed measurement errors. It was found, how- 
ever, that the white dwarf locus is insufficiently distinct in 
the RPMD obtained using these survey parameters. A slight 
restriction of fj.iim for objects with R > 20.5 does much 
to solve the population discrimination problem. Therefore 
the survey parameters given above, with the exception of 
a Hiim of 60 mas/yr for object with R > 20.5, appears to 
be an acceptable compromise between sample maximisation 
and potential contamination and population discrimination 
problems. 

The RPM can be expressed in terms of tangential ve- 
locity ,Vt ,and absolute magnitude, M : 

ii" = M-H51ogioVT- 3.379. (11) 

Evans (1992) has produced theoretical RPMDs by using 
expected 5 logjQ Vt distributions and absolute magnitude 
- colour relations for various populations. Although these 
theoretical predictions were made for specific fields and in- 
corporated error estimates peculiar to that work, they serve 
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Figure 11. RPMD of the survey sample showing the white dwarf 
population divider used. 

as a useful guide to the expected distribution of stellar pop- 
ulations on the RPMD. The RPMD white dwarf population 
locus is found to be an unambiguous population discrimi- 
nator for all colours blue-wards of (O — E) ^ 1.8, with an 
increasing chance of contamination from the spheroid main 
sequence population red-wards of this colour. Transform- 
ing from (O - R) to (B - R) (Evans, 1989 - equation 13), 
spheroid population contamination should become a prob- 
lem red-wards of {B~R) ~ 1.6 which in turn corresponds to 
{Bj — R) ~ 1.2. Indeed on inspection of Figure ^ the white 
dwarf locus does begin to become confused at this colour. 
This contamination occurs solely from the direction of small 
RPM, and the RPMD can still be used with some confidence 
as a population discriminator for objects with high measured 
RPM red-wards of (Bj — R) ^ 1.2. In order to accommodate 
these contamination considerations, the preliminary survey 
white dwarf sample is defined as those objects blue- wards (in 
both plots) of the lines in Figure The lines cut the top of 
the white dwarf locus at (_Bj — i?) ~ 1.2 but allow slightly 
redder objects with higher RPMs into the sample. The white 
dwarf locus is unambiguous blue-wards of (Bj — J?) ~ 1.2. 

Every object in the sample must appear in at least 15 
stacks in each passband. The positional data as a function 
of time have been scrutinised for every object selected as 
a white dwarf candidate, and those with dubious motions 
rejected. While such a process may seem rather arbitrary, 
it was necessary to incorporate this screening stage in the 
sample extraction because simple automated rejection algo- 
rithms such as the 3a rejection routine used here cannot be 
guaranteed to eliminate spurious motions. Some examples 
are shown in Figures |l^ and |l^. All three objects shown 
successfully satisfied all the survey criteria. The object plot- 
ted at the top of Figures |l| and |l| (KX27) shows a clear, 
genuine motion in both x and y in both passbands and was 
included in the final sample without hesitation. The mid- 
dle object (KX18) has larger positional uncertainties and a 
smaller overall motion, but still shows consistent, smooth 
motions and was also included. The final object shows evi- 
dence of large non-linear deviations in the last four epochs of 
the X measures in both passbands. Although the bad-point 
rejection algorithm has removed at least one datum from 
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Epoch Epoch Epoch Epoch 



Figure 12. Some example proper motion plots from the R plate 
stacks. The top plot shows data for object KX27, the middle for 
KX18 and the lower for a rejected object 



each X plot (as shown by the multiple straight line fits), this 
object shows no evidence of proper motion based on the first 
16 data points and certainly cannot be considered a reliable 
proper motion object candidate. This object, along with 9 
others, were rejected from the final WD sample. These re- 
jected objects tended either to have large offsets from a po- 
sitional distribution otherwise consistent with zero motion 
at either the first or last few epochs, as is the case with the 
rejected object described above; or the positional measures 
had an unusually large scatter around the mean position, 
indicating the error in position was larger than the objects 
magnitude would suggest. The reasons for these larger errors 
may be unusual image morphology or the effect of proximity 
to a neighbouring object. While this final rejection proce- 
dure is unsatisfactory in terms of its lack of objectivity, it is 
certainly preferable to inclusion of such objects in the final 
sample, or the introduction of extremely stringent survey 
limits which would doubtless exclude genuinely interesting 
objects. The digitised images have also been inspected. 

The final sample consists of 56 objects which fully sat- 
isfy the photometric, proper motion and RPM/colour survey 
limits. A further two objects, which satisfy the photomet- 
ric and proper motion limits but fall marginally outside the 
RPM / colour cut shown in Figure have also been included 
after favourable follow up observations detailed in Section bI 



Figure 13. Proper motion plots complementary to Figure |l^ 
from the Bj plate stack data 

4 VERY HIGH PROPER MOTION, FAINT 
OBJECT SENSITIVITY LIMITS 

In the previous Section, we discussed in some detail the 
checks made to establish a clean astrometric and photomet- 
ric catalogue. The availability of plates over such a wide 
epoch range as detailed in Table |l| allowed us to search for 
faint and/or very high proper motion objects in this field. 
This is important for a number of reasons. For example, it is 
crucial to firmly establish what the upper limit of detectable 
proper motion is for the methods used, and also to check if 
significant numbers of objects have been missed because of 
this limit or because of the pairing algorithm used. Also, 
it is important to check at fainter magnitudes for very dim, 
high proper motion objects since it is the coolest (and there- 
fore faintest) objects that constrain the age determination 
based on the turn-over in the WDLF. It is also interesting 
to search for very faint, high proper motion halo WDs in the 
light of the current debate concerning the origin of the dark 
lensing bodies detected in microlensing experiments (eg. Is- 
ern et al. 1998 and references therein). 

We performed three experiments to investigate high 
proper motion and/or faint objects: 

(i) Within the restricted epoch range 1992 to 1996, 
ie. five separate epochs, each consisting of a stack of four 
R plates, we used completely independent software employ- 
ing a 'multiple-pass' pairing technique aimed specifically at 
detecting high proper motion objects. This software has suc- 
cessfully detected a very cool, high proper motion degenerate 
WD 0346-1-246 elsewhere (Hambly, Smartt & Hodgkin 1997). 
The pairing algorithm described previously used a 200/xm 
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search radius over 19 yr resulting in an upper limit of 
~ 1 arcsec yr~^, whereas in the multiple-pass test we used a 
maximum search radius of 650/im over a 4 yr baseline, the- 
oretically enabling detection of objects with annual motions 
as high as ~ 10 arcsec. The highest proper motion object 
detected in the catalogue was relatively bright (R ~ 14), 
with fj, ~ 0.8 arcsec yr~^. This experiment revealed two 
/X 0.8 arcsec yr~^ objects , the one mentioned above and 
another slightly fainter object (R ~ 15) ; no objects were 
found with motions larger than this. The colours and re- 
duced proper motions of the two objects indicate that they 
are M-type dwarfs. The second object was undetected in the 
catalogue due to a spurious pairing, an increasingly likely 
scenario for high proper motion objects detected without a 
multiple pass algorithm since they move substantially from 
their master frame position. We note that all the expected 
objects having /i > 0.2 arcsec yr~^ detected in the catalogue 
were also found by this procedure. 

(ii) Using the procedure described in the previous Sec- 
tion, but with a relaxed minimum number of epochs, high 
proper motion objects were searched for. With a 200/im pair- 
ing requirement over a maximum epoch separation of 7 yr, 
the upper limit of proper motion was ~ 1.9 arcsec yr~^. The 
two objects having /i ~ 0.8 arcsec yr~^ found in the previ- 
ous experiment were also recovered here; again, no objects 
were found with motions larger than this. 

(iii) To investigate the possibility of fainter objects, we 
stacked up the R band material in groups of 16 plates at 
epochs 1980, 1983 to 1986, 1987 to 1991 and 1992 to 1996. 
Obviously, over any individual four year period an object 
having a proper motion greater than ~ 1 arcsec yr~^ will 
have an extended image and will not be detected to the same 
level of faintness as a stationary star; nonetheless, the ~ 
0.75™ increase in depth afforded by going from 4 to 16 plate 
stacks (eg. Knox et al. 1998) at least allows us to investigate 
the possible existence of objects having /i ~ 0.5 arcsec yr~^ 
down to R ~ 23 (100% complete to R ~ 22) over an area 
of 25 square degrees. In this experiment, all the objects ex- 
pected from the catalogue were recovered; in addition, one 
star was found having R ~ 20, fj, — 0.47 arcsec yr~^ at 
PA = 179° and RA,DEC = 21h30m8.553s, -44°46'24.09 
(J2000.0). This object is the M-type dwarf 'M20' discov- 
ered in the photometric survey of Hawkins & Bessell (1988) 
and has Bj ~ 23. The faintness in the blue passband is the 
reason that this object is absent from the catalogue. Once 
more, no other high proper motion, fainter stars were found. 

These three experiments allow us to be confident that there 
is no large population of objects having /i <; 1 arcsec yr~^ 
down to faintness limits of R ~ 22 and Bj ~ 23. Further- 
more, the cut-off in the WD sequence seen in the reduced 
proper motion diagrams is real, and not an artefact of in- 
completeness. 



5 FOLLOW UP OBSERVATIONS 

While the RPMD technique is a powerful population dis- 
criminator, it is desirable to obtain follow up observations 
of a sub-set of sample members. The principle motivation 
for this is to explicitly demonstrate the applicability of our 
survey technique by confirming the WD status of the sam- 
ple objects via spectroscopy. Spectroscopic observations red- 
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Figure 14. RPMD of the survey sample with spectroscopically 
observed objects numbered and highlighted. Diamonds denote ob- 
jects identified as subdwarfs or M dwarfs, crosses an immediate 
firm WD identification. Dotted circles represent objects with am- 
biguous spectra (see text). The labelled objects have their spectra 
discussed in detail below. 



wards of the WD cut in the RPMD may also be used to 
investigate the possibility of ultra-cool white dwarfs exist- 
ing in the sample; and as a corollary to this, such observa- 
tions allow clearer population delineation in the RPMD. In 
addition, photometric observations through standard filters 
ensure confidence in photographic-to-standard photometry 
transformations, provide useful independent checks of stel- 
lar parameters (eg T^s) derived from fits to photometry and 
may in future allow alternative population discrimination 
via eg. colour-colour plots. 



5.1 Spectroscopy 

This project was allocated 3 nights of observing time in 1996 
between the 8th and 10th of August, and a further 3 nights 
in 1997 between the 5th and 7th of August on the 3.9m 
Anglo- Australian Telescope. A spectral coverage of 4000- 
7500A was obtained with the RGO spectrograph and 300B 
grating in conjunction with the Tek CCD. Various standards 
were observed throughout each usable night and CuAr lamp 
exposures used for wavelength calibration. The data were 
reduced following standard procedures within the IRAF ^ 
environment. 

The strategy behind these observations was to define as 
clearly as possible the constituents of the lower portion of the 
RPMD, from the blue end of the WD locus down to the ex- 
tremely high H objects below the M stars. The spectroscop- 
ically observed proper motion sample objects are shown in 



IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which is operated by the Association of Universities 
for Research in Astronomy Inc., under contract with the National 
Science Foundation of the United States of America. 
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Figure |l4|. Both the unambiguous bluer region of the WD lo- 
cus and the portion of RPMD lying red-ward of our WD cut 
generally consist of fairly bright objects. These regions have 
been probed spectroscopically, albeit more sparsely than the 
cool WD region consisting of mostly very faint objects. A 
combination of poor weather conditions at the AAT and the 
faintness of our CWD sample has rendered high signal-to- 
noise spectra of these objects unobtainable thus far. 

The bluest objects have clearly defined Hydrogen 
Balmer lines with equivalent widths equal to or in excess of 
those typical for DA white dwarfs of similar colour (Green- 
stein and Liebert, 1990). The objects red- wards of our sam- 
ple cut below the main sequence (MS) show spectra clearly 
distinct from cool WDs. Both subdwarfs and high velocity 
M dwarfs have been identified in this region of the RPMD. 

A star lying near the cut-off region of the RPMD clearly 
showing the absence of strong metal features that would be 
present even in a low-metallicity subdwarf is very likely a 
CWD. We apply this line of argument to our data by select- 
ing from the recent models of Hauschildt et al. (1998) a spec- 
trum appropriate to a subdwarf (metallicity [M/H] > —2) 
of similar effective temperature to a given CWD sample 
object spectrum. The model spectrum is smoothed to the 
approximate resolution of the AAT spectra and multiplied 
through by a synthetic noise spectrum commensurate with 
the AAT CWD spectrum in question. The only features 
likely to be visible after this procedure are ~ 5200A Mgl 
and the ~ 4300A CH G-band, and it is in these regions that 
we look for any evidence that our CWD sample objects are 
in fact subdwarfs. 

This procedure has been undertaken for all CWD candi- 
date spectra. Of these, the spectra of the three objects lying 
within the region of serious potential subdwarf contamina- 
tion are displayed, and we restrict comment on the investi- 
gation of the bluer objects to the statement that none show 
any evidence of subdwarf like spectral features. KX58, dis- 
played in Figure ^ is a convincing CWD candidate, show- 
ing a smooth continuum spectrum with no suggestion of the 
metal features apparent in the model subdwarf spectra. The 
noisier spectrum of KX57 in Figure |l^ also shows no evi- 
dence of subdwarf features, although poorer signal-to-noise 
makes the identification less certain. The object 18d lies sig- 
nificantly red- ward of our CWD cut-off on the RPMD, and 
unfortunately its spectrum (Figure is extremely noisy. 
While it is difficult to draw any conclusions from such poor 
data, the dip at ~ 5200A is a reasonable indication that this 
object is a subdwarf or MS star. Thus the position of 18d on 
the RPMD in conjunction with the spectral data evidence 
means this object does not warrant inclusion in our CWD 
sample. 

To summarise the findings of our spectroscopic survey, 
the only objects showing notable deviation from expected 
WD spectra are the objects below the M dwarf portion of 
the RPMD (diamonds in Figure ^4|) and the more ambiguous 
case of object 18d discussed above. 

5.2 Photometry 

CCD photometry of a subsample of our CWD sample was 
obtained between 29th of July and the 4th of August 1997 
on the Im telescope of the South African Astronomical 
Observatory in Sutherland. Johnson- Cousins V, R, I pho- 
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Figure 15. Spectrum of KX58 with comparison model subdwarf 
spectra 



tometry was obtained for all program stars on the Tek 
(512x512) CCD, with B measures also acquired for suf- 
ficiently bright objects. E-region standards were observed 
continuously through each usable night. Observed magni- 
tudes with associated errors are displayed in Table ^. 

These observed magnitudes provide an independent 
check on the accuracy of the SuperCOSMOS photographic 
photometry, and we use the deviations of mphotographic from 
mccD to obtain errors on the B, V, R, I photographic pho- 
tometry of 0.17, 0.14, 0.13 and 0.16 respectively. The CCD 
photometry also allowed tighter estimates of effective tem- 
perature to be derived for observed objects, although they 
did not provide the hoped for useful constraints on log g. 



6 SAMPLE ANALYSIS 

Previous studies of CWD samples (eg LDM, BRL) have of- 
ten benefited from a comprehensive and wide ranging ob- 
servational database, including high quality spectra, optical 
and IR photometry and parallaxes. These observations, in 
conjunction with detailed WD models, allow determinations 
of stellar parameters such as effective temperature, log g, 
atmospheric composition, mass and bolometric luminosity. 
However, since this is a relatively new project and is con- 
cerned with stars of unusually faint apparent magnitude, 
such a database does not yet exist for this sample. It is 
therefore necessary for us to restrict our analysis, in the 
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Figure 16. Spectrum of KX57 with comparison model subdwarf 
spectra 

first place by exploiting the homogeneity of WD masses by 
assuming a common typical log g for our entire CWD sam- 
ple (the 60 stars with measured log g in BRL have a mean 
surface gravity log g — 8.099±0.044), and secondly by treat- 
ing the atmospheric constituent of each star as an unknown 
parameter whose influence on the resulting WDLF must be 
determined later. 

Bergeron et al. (1995) have published a detailed grid of 
model predictions for Johnson-Cousins U, B, V, R, I (and 
IR) photometry and bolometric corrections as a function 
of effective temperature and log g. Making our assumption 
that log g is always equal to 8, values for effective tempera- 
ture and bolometric luminosity assuming both a H and He 
atmosphere can be calculated for every sample object. 

Fitting for Tcfi is achieved by interpolating the model 
grid at lOK intervals for the colour indices (17 — 5), {B — R), 
(V — R) and {V — I). We then evaluate a^t each Tcs inter- 
val using all available colour indices for the object in ques- 
tion. While the photographic photometry is effective in ade- 
quately constraining Tcff, smaller errors are obtainable with 
the SAAO CCD photometry, and it is used where available. 
The resulting x^, Toff distribution yields a fitted value for 
Tefi and an estimate of the associated error, which can be 
used to read off interpolated values of absolute V magnitude 
and bolometric luminosity. This procedure is performed for 
each object for both an assumed H and He atmosphere. 

Distance moduli obtained from these fits allow calcula- 
tion of tangential velocities via the SuperCOSMOS proper 



5000 6000 7000 
Wavelength(A) 

Figure 17. Spectrum of 18d with comparison model subdwarf 
spectra 

motion measures. The distribution of derived tangential ve- 
locities is consistent with expectations for a sample of Disc 
stars (Evans 1992), and shows no evidence of contamination 
from high velocity halo objects. A summary of the results of 
the fitting procedure is shown in Table ^ including Toh and 
A/boi with associated errors and Kan values. Space density 
values for each object are also presented in Table ^; these 
are discussed in the following section. 

An independent check on the validity of our initial 
CWD sample parameters (SuperCOSMOS photometry and 
astrometry) can be made by comparing known CWD sam- 
ples with our data on the RPMD. Since the RPMD is the 
original means of population discrimination it is also inter- 
esting to overplot known subdwarf samples on our RPMD 
to further address the question of potential contamination. 
BRL published observations of a sample of 110 CWDs, in- 
cluding most of the coolest known degenerates. Extensive 
lists of extreme subdwarfs are less easily obtainable. Ryan 
(1989) used a RPM criterion to extract over 1000 subdwarf 
candidates from the NLTT catalogue. Accurate B and R 
photometry was published for these objects, providing a use- 
ful means of delineating the bluer portion of the WD RPMD 
locus. Monet et al. (1992) identified a subset of 17 extreme 
subdwarfs from their CCD parallax program also involving 
Luyten catalogue stars. Although only V, I photometry is 
published for these stars, we use the photometry published 
in Ryan (1989) to define colour transformations allowing 
the 17 subdwarfs from Monet et al. (1992) to be plotted on 
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Table 4. Derived stellcir parameters for CWD sample assuming either a H or He atmosphere. 
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Object 


B 


V 


R 


I 




CTB 


o-v 




o-l 


KX14(1) 


- 


21.43 


20.78 


20.17 






0.04 


0.09 


0.07 


KX14(2) 


- 


21.47 


20.75 


20.20 






0.04 


0.07 


0.07 


KX15 


19.10 


18.53 


18.08 


17.69 




0.03 


0.02 


0.03 


0.03 


KX20 


16.14 


16.12 


16.16 


16.22 




0.02 


0.02 


0.02 


0.02 


KX22 


- 


19.94 


19.66 


19.55 






0.03 


0.04 


0.04 


KX23 


19.52 


19.05 


18.74 


18.61 




0.03 


0.01 


0.01 


0.03 


KX27 


- 


18.15 


17.97 


17.91 






0.02 


0.02 


0.03 


KX29 


18.67 


18.36 


18.17 


17.90 




0.03 


0.02 


0.03 


0.03 


KX30 


- 


21.12 


20.64 


20.19 






0.03 


0.05 


0.08 


KX33 


- 


21.22 


20.69 


20.10 






0.04 


0.06 


0.07 


KX39 


- 


21.19 


20.66 


20.18 






0.04 


0.06 


0.06 


KX41a 


20.26 


19.38 


18.93 


18.41 




0.05 


0.03 


0.04 


0.04 


KX41b 


20.28 


19.33 


18.81 


18.36 




0.04 


0.03 


0.04 


0.04 


KX44 




21.37 


20.81 


20.42 






0.05 


0.07 


0.08 


KX53 




20.96 


20.70 


20.21 






0.03 


0.05 


0.06 


KX57 




21.09 


20.40 


19.84 






0.04 


0.06 


0.08 



Table 3. Johnson-Cousins CCD photometry taken at SAAO for 
selected members of our CWD sample. Object KX41 was resolved 
as a double-degenerate on the CCD frame and thus has photom- 
etry for each component. Object KX14 had two independent sets 
of observations, shown as (1) and (2) 
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Figure 18. A comparison of various populations on the RPMDs: 
crosses denote our CWD sample , circles the BRL CWD sample, 
asterisks the Monet et al. extreme red subdwarfs and dots the 
brighter Ryan subdwarfs 



assess contamination of the CWD sample from the direc- 
tion of small RPM. However the subdwarf RPM locus is not 
predicted to deviate significantly from a straight line in the 
CWD colour regime (Evans 1992), and if we take the high H 
extent of the two subdwarf samples plotted to be indicative 
of the limit of the extreme subdwarf locus on the RPMD, 
the dashed lines plotted on Figure |l^ should be a good guide 
to the limit of the subdwarf locus for the intermediate colour 
range. It may then immediately be seen that the vast ma- 
jority of our CWD sample is safely within the WD region of 
the RPMD. The two redder borderline stars have reasonable 
spectroscopic confirmation of their WD status (Figures |l^ 
and |l^), leaving only one potentially dubious object. 



the {B — R), RPMD planes. These objects define a portion 
of the RPMD marginally red-wards of the faintest CWDs 
where the most extreme contaminants may be expected to 
lie. 

Figure ^ shows the two RPMD with the four sam- 
ples plotted. There are several points to be made concern- 
ing this plot. Firstly, our CWD sample and the the BRL 
sample of previously known CWDs lie on the same region of 
the diagram, providing further confirmation of the validity 
of our survey procedure. It can also be seen that the BRL 
sample contains redder, cooler stars than our sample. This 
may be expected since the BRL sample is rather eclectic 
and contains some of the coolest WDs known, whereas our 
sample is drawn from a rigidly defined survey in a partic- 
ular ESO/SERC field. We note also that the cool portion 
of the BRL sample does not extend into the portion of the 
RPMD beyond our population discrimination cut-off shown 
in Figure ^ which may be interpreted as indicating that we 
are not failing to sample portions of the RPMD containing 
CWDs (but see Section M below). Both subdwarf samples 
lie in clearly distinct regions of the RPMD to our sample, 
although the cooler subdwarfs are all too red to directly 



7 THE WDLF 

In order to construct a WDLF, space densities must be cal- 
culated for a survey limited by both apparent magnitude and 
proper motion. The CWD survey sample presented in Ta- 
ble ^consists of stars with widely varying intrinsic brightness 
and tangential velocity, and is therefore not volume-limited 
(since for example intrinsically bright objects are sampled 
out to greater distances than the coolest, faintest stars). The 
standard solution to this problem, Schmidt's (1968, 1975) 
1/Vmax estimator, has been extensively studied with specific 
reference to the WDLF (Wood & Oswalt 1998). The 1/Vmax 
method assigns each sample object a survey volume defined 
by the maximum distance dmax the object could have and 
satisfy the survey limit criteria. For this survey, an object 
at distance d with proper motion /i and magnitudes B and 

R has dmax 



dmax = min 



d- 



,dlO 



-R) 



dm 



,0.2(Sh 



-B) 



(12) 



For the simple case of uniform stellar space density, the sur- 
vey field solid angle Q, can then be used to calculate the 
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corresponding Vmax {~ Q.d^^^/'i) . However, many of the 
objects in our CWD sample yield a dmax comparable to the 
scale height of the Disc, and we follow the method of Tin- 
ney, Reid & Mould (1993) in generalising our calculation of 
Vmax to allow for the truncation of the survey volume by 
the scale height effect. In this prescription, the volume is 
found by integrating over an exponentially decreasing den- 
sity law (Stobie, Ishida & Peacock 1989) with scale height 
h at galactic latitude b, yielding the modified expression for 
the volume out to a distance d: 



V = n-^{2-(^2 + 2C + 2)e-«}, 



(13) 



where ^ = dsin b/h (Equation 9 in Tinney et al. 1993). Each 
object is thought of as a sampling of the survey volume 
Vmax, and thus contributes a space density of V^^ax to its 
particular LF bin. We adopt the convention of LDM in as- 
signing the uncertainty in each space density contribution as 
being equal to that contribution (ie. 1 ± 1 objects per sam- 
pling). We may therefore construct a LF simply by summing 
the individual space density contributions in each luminosity 
bin, and the error obtained by summing the error contribu- 
tions in quadrature. A true reflection of the LFs observa- 
tional uncertainties should also allow for the uncertainties 
inherent in photometric fits to model atmospheres described 
in the previous section. The errors in Mboi detailed in Ta- 
ble 1^ suggest the introduction of horizontal error bars in any 
observational LF is necessary. When converting our sample 
Mboi values into luminosity units via 



Mb. 



-2.5 log L/Lq 4.75 



(14) 



we also calculate the upper and lower 1 sigma luminosity 
uncertainties for each object. For a bin containing N objects 
we then combine eg. the upper luminosity 1 sigma uncer- 
tainties, (T„, using 



N o 



N 



(15) 



to yield an estimate for the horizontal error bar au, with an 
analogous procedure for the lower luminosity error bounds. 
In addition, to give the most realistic estimate for the LF in 
magnitude bins containing a very few objects, we plot the 
binned data at the mean luminosity of the objects in the 
bin, rather than at the mid-point of that bin. 

Table ^ gives the LF calculated in this fashion for in- 
teger magnitude bins. Only the cool end of the LF is given 
here (Mboi > 12.25), and a field size of 0.0086 steradians 
and Disc scale height of 300 parsecs were used in calculating 
space densities. Columns 3 and 4 give the plotted (Figure |l9[) 
LF with upper and lower error bounds inparenthesis. The 
hot WD data point not detailed in Table n represents the 6 
stars with (3 > — log L/Lq > 2). These hotter objects tend 
to have large errors in fitted Mboi, and the resulting large 
horizontal error bars makes broader binning appropriate. It 
is necessary to choose either a pure hydrogen or helium at- 
mosphere for each object to construct the LF. We use our 
'best guess' atmospheres for the LF described here: for each 
object we choose the atmosphere with the lower model 
fit to the photometry, with the exception of objects with 
6000 > Tcffn > 5000 which are deemed to be occupying 
the 'non-DA gap' at this temperature (BRL) and are there- 
fore automatically designated a pure H atmosphere. Note 



CO 



o 

m 
o 



-Log(L/LQ) 

Figure 19. Observational LF with comparisons: this work (filled 
circles), hot WDLF based on Fleming, Liebert and Green (1986) 
from LDM (open circles) and the LRB redetermination of the 
LDM CWDLF (open triangles) 



however that the photometry does not adequately constrain 
atmospheric composition and this 'best guess' LF is only 
one arbitrary realisation of the data (see Section ^ below) . 

The standard means of assessing the completeness of a 
sample analysed using the 1/Vmax method is to calculate 
the mean value of the ratio of Vots, the volume out to an 
object, to Vmax- A complete survey evenly sampling the sur- 
vey volume should yield (Vobs/Vmax) = |. This is generally 
not the case for published CWD samples, although some au- 
thors have incorporated completeness corrections into their 
analyses to account for the effects of the original survey in- 
completeness (OSWH). The (Vots/Vmax) calculated for this 
sample is 0.495 or 0.496 choosing either all H or He atmo- 
spheres. From a complete sample containing 58 objects we 
expect (Vobs/Vmax) = 0.5 ± 0.038, indicating our sample 
is consistent with being drawn from a complete survey. It 
should be emphasised however that this result cannot be re- 
garded as proof of completeness, since clearly an incomplete 
survey sample may also exhibit (Vobs/Vmax) ~ ^■ 

The total space density determined from our 'best 
guess' sample is 4.16 x 10~^ WDs per cubic parsec, approx- 
imately 25% greater than that found by LRB. These results 
are certainly consistent, since the simulations of Wood and 
Oswalt (1998) predict errors of ~ 50% in total space den- 
sity estimates from samples of 50 CWDs using the 1/Vmax 
technique and additional uncertainties are introduced by our 
lack of knowledge of the WD atmospheric constituents. Our 
findings reiterate that WDs represent only a small fraction 
(~ 1%) of the local dynamically estimated mass density. In- 
terestingly, we do not confirm the much higher total WD 
space densities found recently by two independent studies. 
We note however that these studies (Ruiz & Takamiya 1995, 
Festin 1998) make only a tentative claim to detection of a 
high WD space density due to the small samples (N < 10) in- 
volved. A third study (OSWH) searched exclusively for WDs 
in CPMBs, and found a total space density of 5.3 x 10^'^ for 
these objects. At this space density, and given our survey 
area and the resolution of the COSMOS data, we would not 
expect to find any CPMB WDs. It is therefore not surpris- 
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Bin Center 




Space Density 




Mboi - log L/Lq 


-logL/Lg 




Number 


12.75 3.20 


3.27(3.05,3.45) 


-3.39(-3.23,-3.65) 


14 


13.75 3.60 


3.62(3.44,3.76) 


-3.32(-3.18,-3.51) 


13 


14.75 4.00 


4.04(3.92,4.15) 


-2.68(-2.56,-2.84) 


20 


15.75 4.40 


4.30(4.18,4.40) 


-2.97(-2.80,-3.24) 


5 



Table 5. Luminosity Function 



ing that no object in our sample is a CPMB member, since 
the survey technique is only sensitive to lone WDs or double 
degenerate binary systems. 



8 DISCUSSION 

An estimate of the Disc age may be obtained by comparing 
our WDLF with expectations from theoretical models. We 
compare our data with two sets of models, which are avail- 
able in the form of curves at integer 1 Gyr Disc age intervals. 
Since we are fitting to a cut-off in space density, the lack of 
detected objects beyond our faintest observational bin as- 
sumes added significance. We can calculate the probability 
of detecting zero objects in the next faintest luminosity bin 
because we have well defined survey limits: at given faint 
luminosity the proper motion survey limit is irrelevant and 
the survey is sampling a known volume defined by the pho- 
tometric survey limits. This volume was calculated using B 
and R magnitudes for very cool WDs from the recent mod- 
els of Ifansen (1998), which combined with the photometric 
survey limits yield a dmax for both an If and an He atmo- 
sphere WD. The minimum dmax defines the survey volume 
at that magnitude. We have assumed the LF at any given 
luminosity consists of equal numbers of H and He WDs, in 
which case the H WD survey volume provides the best con- 
straint on the LF and is adopted for fitting. Using Poisson 
statistics, the probability of detecting zero objects in a vol- 
ume V in which a model predicts space density p is simply 
e"''^. We assume the errors on the data points are normally 
distributed, and derive best fits by maximising the proba- 
bility of the dataset for each curve and comparing the fits 
for the various Disc ages. 

The various inputs to the first set of WDLFs we use, 
those of Wood, were described in detail in Wood (1992). 
Careful consideration was given to the various inputs to 
the WDLF, such as initial mass function, star formation 
rates and initial-final mass relation. The Wood WD evo- 
lutionary models for the WDLFs have since been updated 
(Wood 1995, OSWH), and consist of mixed carbon-oxygen 
core WDs with hydrogen and helium surface layer masses 
of 10^* and lO^^M* respectively, and also utilise revised 
opacities and neutrino rates. It is these more recent model 
WDLFs that are used here. 

The best fit of our 'best guess' LF to the Wood model 
WDLFs is for a Disc age of 9 Gyr, and is shown in Figure po[ 
Although this gives a reasonable first indication of the Disc 
age implied by our sample, further investigation is neces- 
sary since the photographic photometry fitting procedure 
used to estimate Tea (as described in Section |^ does not re- 
liably constrain the atmospheric constituent of the stars in 
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Figure 20. The Wood WDLF with a Disc age of 9 Gyr - the 
best fit to our 'best guess' LF. 8 Gyr (below) and 10 Gyr (above) 
curves are also shown in dashed lines. 



our sample. We have addressed this question by construct- 
ing a large number of WDLF 'realisations' from our sample 
data, each time giving each star a 50% probability of having 
either a H or He atmosphere. Every resulting WDLF was fit 
to the models and the best fit Disc age recorded. We have 
also used this analysis to assess the effect of binning on the 
fitting procedure. Figure ^ (a) displays the results of fits to 
1000 realisations binned in 1.0 Mboi bins and a further 1000 
in 0.75 Mboi bins. These results give a fairer indication of 
the Disc age and associated errors than Figure which is 
effectively one arbitrarily picked realisation. 

The second set of theoretical WDLFs we use are de- 
scribed in Garci'a-Berro et al. (1997) (henceforth GB mod- 
els). These LFs include the expectation that the progeni- 
tors of the faintest WDs are likely to have been massive 
stars since these stars evolve more quickly and the resulting 
(Oxygen-Neon) massive WDs also cool faster. These models 
also include a predicted delay in Carbon- Oxygen WD cool- 
ing induced by the separation of C and O at crystallization 
(Hernanz et al. 1994). The incorporation of these considera- 
tions into the theoretical WDLF leads to a broader predicted 
peak at the tail of the LF. We find a best fit of 11 Gyr to 
our sample, as shown in Figure Again, we have investi- 
gated the effect of our poor knowledge of our samples at- 
mospheres in the same way as above. The results, shown in 
Figure ^ (b), highlight the effect that variations in binning 
can have on the fits for this second set of models. 

The major contributor to this problem seems to be that 
the model curves are sufficiently indistinct in the region of 
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Figure 21. Results of constructing WDLFs from our sample by 
randomly assigning H or He atmospheres. Full line histograms 
denote 1 M^oi bins, dotted lines 0.75 M^oi bins. For each bin size 
1000 realisations of the WDLF were constructed and fit to the 
models - the number of times each Disc age was the best fit is 
displayed. Graph (a) shows fits to the Wood models, graph (b) 
fits to GB. 




Figure 22. The best fit of our 'best guess' LF to the GB WDLF: 
a Disc age of 11 G yr w ith 'adjacent' Disc age LFs shown in dashed 
lines as in Figure koI. 

our sample dataset (as can be seen clearly in Figure p^ ) that 
changes in observational LF binning can have a significant 
bearin g o n the result of the fit. This effect is apparent in 
Figure Ell (b), where a bin size of 0.75 Mboi yields a strong 
preference of a Disc age of 14 Gyr, in contrast to the 12- 
13 Gyr Disc age selected by the other binning regime (note 
also that 14 Gyr was the oldest curve available for fitting). 
Although it is difficult to resolve this matter satisfactorily 
with the current data set, there are two pertinent points to 
raise. The problem does not apply to the Wood models, the 
fits to which are extremely difficult to dislodge from the 9- 
11 Gyr region indicated by Figure |l| (a). Secondly, Wood 
and Oswalt (1998), as a result of their Monte Carlo analysis 
of the WDLF, recommend choosing a binning in which the 
crucial lowest luminosity bin contains ~ 5 objects; giving 



a reasonable compromise between the requirements of good 
signal to noise in the final bin and having that bin as faint 
as possible to provide maximum information on the position 
of the cut off. This means a bin size of 1 Afboi for our data 
set, for which the fitted Disc age is well constrained for both 
models. 

A brief investigation of the effect of altering the Disc 
scale height in Equation |l3| revealed that for any scale height 
between 250 and 400 pc the alterations in space densities for 
a LF binned as in Table ^ are restricted to a few hundredths 
in log (f). Variations at this level have a negligible effect on 
the fits to model LFs. This behaviour is to be expected , 
since it is objects with large dmax which are most affected 
by changes in the scale height, and these make the smallest 
contributions to the space density when using the 1/Vmax 
prescription. 

An overview of all the various fittings to both sets of 
models points firmly to a Disc age between 9 and 11 Gyr us- 
ing the Wood models, and an age of 12-14 Gyr using the GB 
models. This discrepancy is expected (LRB, Garcia-Berro 
et al. 1997) and indicates the extent of errors introduced by 
uncertainties in the WD crystallisation process; the Wood 
model Disc age should be a reliable lower limit however, and 
Figure |2l| demonstrates the difficulty in obtaining a Disc age 
of less than 9 Gyr from our data. We note that the Wood 
models seem to represent our data better, the GB models 
not following the peak in the observed LF. Wood (1992) re- 
ported that the uncertainties in the inputs to model WDLFs 
lead to a further ~ 1 Gyr contribution to the total error. Fur- 
ther insight into the error associated with the Disc age may 
be gained by considering the extremities of the distribution 
of atmosphere types within the sample. The 50% probabil- 
ity assigned to the H and He atmosphere is adopted in the 
absence of strong evidence for a dominant atmospheric type 
amongst the crucial coolest WDs (see eg. the coolest bin in 
Figure 1 of BRL), and has no physical basis. Adopting all 
H atmospheres leads to high space density estimates in the 
final LF bins and arbitrarily old Disc age estimates. More 
interestingly, a LF composed of all He atmosphere WDs is 
still incompatable with Disc ages below 8 Gyr, regardless of 
the model WDLF used. This reiterates the important point 
that the cosmologically interesting lower Disc age limit ap- 
pears to be 8 Gyr, and that even a Disc as young as this 
must be considered unlikely (Figure ^]). Our adopted Disc 
age estimation is therefore IOI1I3 Glyr. 

There are at least two possible effects arising from our 
present lack of extensive follow up observational data that 
could affect the derived Disc age. First, as may be seen in 
Figure ^ there is a small group of objects lying just beyond 
our RPMD cut-off that have not been included in our CWD 
sample. Although the distribution of known WDs on the 
RPMD (the BRL objects in Figure |l|) indicates a popula- 
tion of CWDs is not expected in this region, the reason for 
this may be that many of the known CWDs have been them- 
selves selected on the basis of RPM criteria (Luyten 1970). 
It may be that the region of the RPMD just beyond our 
cut-off has not been adequately investigated before, since in 
the presence of noisier data from blink comparators and eye- 
measure photometry it would be hopelessly confused. These 
considerations argue strongly for a wholesale spectroscopic 
survey using multi-fibre instruments of the entire population 
below Hr = 19 for complete confidence that our sample does 
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not exclude any CWDs. For the present, the possibility that 
a few more CWDs exist in our catalogue but do not sat- 
isfy our RPMD survey criteria cannot be ruled out. Such 
objects would certainly be very cool, resulting in a higher 
Disc age estimate. The second effect concerns the question 
of mass. The difficulty is amply demonstrated by the case of 
CWD ESQ 439-26 (Ruiz et al. 1995), which was observed to 
have a luminosity fainter by 1 magnitude than the WDLF 
cut off. Analysis of the object's optical energy distribution 
in conjunction with a measure of trigonometric parallax al- 
lowed the authors to conclude that it's low luminosity was 
in fact due to its large mass, or small radius. Again, ideally 
it would be desirable to obtain parallaxes and CCD opti- 
cal and IR photometry of our entire sample, obviating the 
need to assume a mass of 0.6Mq (and allowing the surface 
composition to be constrained). 

In summary, our WD sample has passed every test for 
completeness applied to it. The calculated space density of 
WDs is slightly higher than that for the LDM sample, how- 
ever we do not detect the much higher total space densities 
found by more recent authors (Ruiz 1995, Festin 1998). Our 
WDLF yields a Disc age estimate of lO^g Gyr, older but 
still consistent with previous estimates (Winget 1987, LRB, 
OSWH). In the context of current cosmochronometry, our 
Disc age estimate is consistent with current Globular clus- 
ter age estimates of 13-14Gyr (Vandenberg 1998). Finally, 
when combined only with a conservative 1 Gyr value for 
the halo-Disc formation interval (Burkert, Truran & Hensler 
1992, Pitts & Tayler 1992) and a further 1 Gyr for the big 
bang-halo formation interval, a 10 Gyr Disc excludes = 1, 
A = cosmologies based on current estimates for Ho of GO- 
SO (Freedman 1998). 



9 CONCLUSIONS 

(i) Using a large collection of COSMOS/SuperCOSMOS 
digitised Schmidt plate data in ESO/SERC field 287 we have 
extracted a sample of proper motion objects. Number counts 
indicate this sample is complete down to the survey proper 
motion limit, which was chosen with care to exclude con- 
taminant spurious proper motion measures. 

(ii) A sample of cool white dwarfs have been culled from 
our proper motion objects using the reduced proper mo- 
tion technique. By overplotting samples of known extreme 
subdwarfs we show our CWD sample is unlikely to be con- 
taminated by other stellar population groups. We have con- 
firmed the WD status of a number of our sample with AAT 
spectroscopy. The sample passes the (V/Vmax) completeness 
test. 

(iii) We calculate a total WD space density of 4.16 x 10~^ 
WDs per cubic parsec using Schmidt's (l/Vmax) method. 
Careful comparison of luminosity functions constructed from 
our sample and theoretical models indicate an age for the lo- 
cal Galactic Disc of lO^g Gyr, older than previous estimates 
using this technique. 
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APPENDIX A: REDUCED PROPER MOTION 
DIAGRAMS 

The reduced proper motion (RPM) is defined by 

// = m + 51ogio^i + 5 (Al) 

where m is apparent magnitude and ^ proper motion. A re- 
duced proper motion diagram (RPMD) is a plot of colour 
against RPM. It is a powerful way of combining proper 
motions and photometry to distinguish stellar population 



groups. Equation Al can be re-written using the relation- 



ships m = M — 5 + 51ogd and /i = VT/4.74ci to give 

H = M + ^ logio Vt - 3.379 (A2) 

where M is the absolute magnitude, Vt the transverse ve- 
locity (in kms~^) and d the distance. Since M and Vt are 
both intrinsic properties of the star, so too is H. 

To see the significance of H, suppose that every star 
had an identical Vt; H would then clearly be simply M plus 
a constant and the distribution of a particular population 
group in H at a particular colour would depend solely on 
the spread of the populations' colour-magnitude relation at 
that colour. Of course there is a distribution in 51ogj^Q Vt 
for each population, but since the tangential velocities are 
distributed around a most probable value the RPM serves as 
an estimate of M, ie. M = a bH (a and b constants). The 
resulting locus for each population in the RPMD is then 
the convolution of its colour magnitude distribution with 
its 5 logjQ Vt distribution over the diagrams colour range. 
To allow population discrimination in some colour region 
we therefore only require that the various population loci 
do not overlap significantly in that region of the RPMD. 
In effect the RPMD is analogous to the Hertsprung-Russell 
diagram, and in both plots the white dwarf population is 
quite distinct in most colours. 
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